arXiv:nucl-th/0211015vl 6 Nov 2002 


Low Mass Neutron Stars and the Equation of State of Dense Matter 


J. Carriere * and C. J. Horowitz^ 

Nuclear Theory Center and Dept, of Physics, Indiana University, Bloomington, IN f7f05 

J. Piekarewicz 1 

Department of Physics, Florida State University, Tallahassee, FL 32306 
(February 4, 2008) 

Neutron-star radii provide useful information on the equation of state of neutron rich matter. 
Particularly interesting is the density dependence of the equation of state (EOS). For example, the 
softening of the EOS at high density, where the pressure rises slower than anticipated, could signal 
a transition to an exotic phase. However, extracting the density dependence of the EOS requires 
measuring the radii of neutron stars for a broad range of masses. A “normal” 1.4Mq (MQ=solar 
mass) neutron star has a central density of a few times nuclear-matter saturation density (po)- In 
contrast, low mass (~ 0.5 Mq) neutron stars have central densities near po so its radius provides 
information on the EOS at low density. Unfortunately, low-mass stars are rare because they may 
be hard to form. Instead, a precision measurement of nuclear radii on atomic nuclei may contain 
similar information. Indeed, we find a strong correlation between the neutron radius of 208 Pb and 
the radius of a 0.5M 0 neutron star. Thus, the radius of a 0.5 Mq neutron star can be inferred from 
a measurement of the the neutron radius of 208 Pb. Comparing this value to the measured radius of 
a ~1.4Mg neutron star should provide the strongest constraint to date on the density dependence 
of the equation of state. 


I. INTRODUCTION 

The structure of neutron stars, particularly their 
masses and radii, depend critically on the equation of 
state (EOS) of dense matter Q]. New measurements of 
masses and radii by state-of-the-art observatories should 
place important constraints on the EOS. Observing a 
rapid change of the EOS with density could signal a tran¬ 
sition to an exotic phase of matter. Possibilities for new 
high density phases include pion or kaon condensates [|| , 
strange quark matter j3) , and/or a color superconduc¬ 
tor Measuring neutron-star radii R(M) for a large 

range of neutron star masses M is attractive as it would 
allow one to directly deduce the EOS Q], that is, the 
pressure as a function of the energy density P(e). 

While the masses of various neutron stars are accu¬ 
rately known j^], precise measurements of their radii do 
not yet exist. Therefore, several groups are devoting con¬ 
siderable effort at measuring neutron-star radii. Often 
one deduces the surface temperature T 0 0 and the lumi¬ 
nosity L of the star from spectral and distance mea¬ 
surements, respectively. Assuming a black-body spec¬ 
trum, these measurements determine the surface area, 
and thus the effective radius Roc, of the star from the 
Stefan-Boltzmann law: 

L = 4tt aRlT^ . (1) 

Opportunities for precision measurements on neutron- 
star radii include the isolated neutron star RX J185635- 
3754 g,|,|T^| and quiescent neutron stars in globular clus¬ 
ters, such as CXOU 132619.7-472910.8 [|lT), where dis¬ 
tances are accurately known. Moreover, Sanwal and col¬ 
laborators have recently detected absorption features in 


the radio-quiet neutron star IE 1207.4-5209 jlj] that may 
provide the mass-to-radius ratio of the star through the 
determination of the gravitational redshift of the spec¬ 
tral lines. These observations are being complemented 
by studies that aim at constraining the composition of 
the neutron-star atmosphere iQ. Finally, models of ro¬ 
tational glitches place a lower limit on the radius of the 
Vela pulsar at Roo g) 12 km [ pT[ . 

While a determination of the mass-radius relation 
R(M) for a variety of neutron stars would uniquely 
determine the equation of state, unfortunately all ac¬ 
curately determined masses to date fall within a very 
small range. Indeed, a recent compilation by Thorsett 
and Chakrabarty of several radio binary pulsars place 
their masses in the narrow range of 1.25 — 1.44 Mq 0. 
Note that several X-ray binaries appear to have larger 
masses, perhaps because of accretion. These include 
Cyg X-2 with a mass of 1.8 ± 0.2M 0 flul , Vela X-l with 
1.9 M 0 pl l, and 4U 1700-37 |TJ]. If confirmed, they 
could provide additional information on the high density 
EOS. However, these mass determinations are not with¬ 
out controversy |l^,|TfJ. On the other hand, it may be 
difficult to form low mass neutron stars from the collapse 
of heavier Chandrasekhar mass objects. If so, informa¬ 
tion on the low density EOS may not be directly available 
from neutron stars. Thus, it is important to make max¬ 
imum use of any R(M) measurements even if these are 
available for only a limited range of masses. 

Additional information on the low density EOS may 
be obtained from precision measurements on atomic nu¬ 
clei. For example, the neutron radius of a heavy nu¬ 
cleus, such as 208 Pb, is closely related to the pressure 
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of neutron rich matter p0| . Indeed, heavy nuclei de¬ 
velop a neutron-rich skin in response to this pressure. 
The higher the pressure the further the neutrons are 
pushed out against surface tension, thereby generating 
a larger neutron radius. However, nuclear properties de¬ 
pend only on the EOS at normal (in the interior) and 
below (in the surface) nuclear-matter saturation density 
(po ss 0.15 nucleons/fm 3 ). This is in contrast to conven¬ 
tional 1.4 Mq neutron stars that, with central densities 
of several times po, also depend on the high-density com¬ 
ponent of the EOS. This is not the case for low mass 
neutron stars (of about 1/2 Mq). Reaching central den¬ 
sities near po, l°w mass neutron stars probe the EOS at 
similar densities as atomic nuclei. Therefore, one could 
infer the radius of low mass neutron stars from detailed 
measurements on atomic nuclei. Having inferred the radii 
of low mass neutron stars from a nuclear measurement, 
combined with the measured radius of a IAMq star, may 
enable one to deduce the density dependence of the EOS. 
(For a recent discussion on the minimum stable mass of a 
neutron star see Ref. .) The parity radius experiment 
at Jefferson Laboratory |22|| aims to measure accurately 
and model independently the root-mean-square neutron 
radius ( R n ) of 208 Pb via parity violating elastic electron 
scattering [|U . Such an experiment probes neutron den¬ 
sities because the weak vector charge of a neutron is much 
larger than that of a proton. The goal of the experiment 
is to measure R n to a 1% accuracy (within ss ±0.05 frn). 

The outline of the paper is as follows. In Sec. [H] we 
present a relativistic effective field theory formalism to 
study relationships between the neutron radius of 208 Pb 
and the radii of neutron stars. Uncertainties in these re¬ 
lationships are estimated by considering a wide range of 
effective field theory parameters, all of them constrained 
by known nuclear properties. This formalism has been 
used previously to study correlations between the neu¬ 
tron radius of 208 Pb and the properties of the neutron- 
star crust HU, the radii of 1.4 Mq neutron stars Hs| , 
and the direct URCA cooling of neutron stars HU. As 
the radius of low mass neutron stars depends on the solid 
crust of nonuniform matter, a treatment of the crust is 
discussed in Sec. HI. Our results, presented in Sec. 0, 
show a strong correlation between the radius of a low 
mass neutron star and the neutron radius R n of 208 Pb 
that is essentially model independent. This is because 
the structure of both objects depend on the EOS at sim¬ 
ilar densities. In contrast, the radius of a 1.4 Mq neu¬ 
tron star shows a considerable model dependence. This 
is because a 1.4 Mq neutron star is also sensitive to the 
EOS at higher densities, and the high density EOS is 
not constrained by nuclear observables. In Sec. ^ we 
conclude that properties of low mass neutron stars can 
be inferred from measuring properties of atomic nuclei. 
In particular, the radius of a 1/2 Mq neutron star can 
be deduced from a measurement of R n in 208 Pb. One 
will then be able to directly compare this inferred ra¬ 


dius to the measured radius of an « 1.4 Mq neutron 
star to gain information on the density dependence of 
the EOS. Thus, even if a mass-radius measurement for a 
single (~ 1.4 Mq) neutron star is available, one can use 
the atomic nucleus to gain information on the density 
dependence of the EOS. This should provide the most 
precise determination of the density dependence of the 
EOS to date and should indicate whether a transition to 
a high density exotic phase of matter is possible or not. 


II. FORMALISM 


Our starting point is the relativistic effective-field the¬ 
ory of Ref. |27| 1 supplemented with new couplings be¬ 
tween the isoscalar and the isovector mesons. This allows 
us to correlate nuclear observables, such as the neutron 
radius of 208 Pb, with neutron star properties. We will 
explore uncertainties in these correlations by consider¬ 
ing a range of model parameters. The model has been 
introduced and discussed in detail in several earlier ref¬ 
erences yet a brief summary is included here 

for completeness. 

The interacting Lagrangian density is given by [|2 l|j27| 


Ant = -0 5 , s0-(s'vV/±yT • b^ + ^(1 +r 3 ) y 


ip 


~ |j(3s0) 3 -^(5s0) 4 ±^.gv(^M 1/M ) 2 

+ b M • b M [A s g 2 0 2 + A v g 2 V/U M ] . (2) 


The model contains an isodoublet nucleon field (0) in¬ 
teracting via the exchange of two isoscalar mesons, the 
scalar sigma (0)and the vector omega (U M ), one isovector 
meson, the rho (b^ 1 ), and the photon ( A M ). In addition 
to meson-nucleon interactions the Lagrangian density in¬ 
cludes scalar and vector self-interactions. Omega-meson 
self-interactions 0 soften the equation of state at high 
density. Finally, the nonlinear couplings A s and A v are 
included to modify the density-dependence of the sym¬ 
metry energy a sym (p) H®!!]- We employ Eq. (|) in 
a mean field approximation where the meson fields are 
replaced by their ground state expectation values. The 
couplings constants in Eq. (j|j) are fit to nuclear mat¬ 
ter and finite nuclei properties. All of the parameter 
sets considered here, namely, NL3 [|U, S271 HU> and 
Z271 HU reproduce the following properties of symmet¬ 
ric nuclear matter: saturation at a Fermi momentum of 
kF = 1-30 fm _1 with a binding energy per nucleon of 
— 16.24 MeV and an incompressibility of K = 271 MeV. 
The various parameter sets differ in their effective masses 
at saturation density, in their w-meson self interactions 
(which are included for Z271 and neglected for NL3 and 
S271) and in the nonlinear couplings A s and A v (see 
Table |l|). Note that the NL3 parametrization has been 
used extensively to reproduce a variety of nuclear prop¬ 
erties @. 
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The symmetry energy at saturation density is not well 
constrained experimentally. However, an average of the 
symmetry energy at saturation density and the surface 
symmetry energy is constrained by the binding energy 
of nuclei. Thus, the following prescription has been 
adopted: the value of the NNp coupling constant is ad¬ 
justed so that all parameter sets have a symmetry energy 
of 25.67 MeV at k F = 1.15 fm -1 . This insures accurate 
binding energies for heavy nuclei, such as 208 Pb. Follow¬ 
ing this prescription the symmetry energy at saturation 
density is predicted to be 37.3, 36.6, and 36.3 MeV for 
parameter sets NL3, S271, and Z271, respectively (for 
A s = A v = 0). Changing A s or A v changes the density 
dependence of the symmetry energy by changing the ef¬ 
fective rho-meson mass. In general increasing either A s 
or A v causes the symmetry energy to grow more slowly 
with density. 

The neutron radius of 208 Pb depends on the density de¬ 
pendence of the symmetry energy. A large pressure for 
neutron matter pushes neutrons out against surface ten¬ 
sion and leads to a large neutron radius. The pressure de¬ 
pends on the derivative of the energy of symmetric matter 
with respect to density (which is approximately known) 
and the derivative of the symmetry energy, da sym /dp. 
Thus parameter sets with a large da sym /dp yield a large 
neutron radius in 208 Pb. Note that all parameter sets 
approximately reproduce the observed proton radius and 
binding energy of 208 Pb. Therefore changing A s or A v 
allows one to change the density dependence of the sym¬ 
metry energy da sym /dp , while keeping many other prop¬ 
erties fixed. Once the model parameters have been fixed, 
it is a simple matter to calculate the EOS for uniform 
matter in beta equilibrium, where the chemical poten¬ 
tials of the neutrons p n , protons p p , electrons p e , and 
muons p p satisfy, 

pn Pp = Pe = Pp • (3) 

Note that the high density interior of a neutron star is 
assumed to be a uniform liquid; possible transitions to a 
quark- or meson-condensate phase are neglected. 

TABLE I. Model parameters used in the calculations. The 
parameter k and the scalar mass m s are given in MeV. The 
nucleon, rho, and omega masses are kept fixed at M = 939, 
m p = 763, and = 783 MeV, respectively — except in the 
case of the NL3 model where it is fixed at m^ = 782.5 MeV. 


Model 

m s 

9s 

gl 

K 

A 

C 

NL3 

508.194 

104.3871 

165.5854 

3.8599 

-0.0159049 

0 

S271 

505 

81.1071 

116.7655 

6.68344 

-0.01580 

0 

Z271 

465 

49.4401 

70.6689 

6.16960 

+0.156341 

0.06 


III. BOUNDARY BETWEEN CRUST AND 
INTERIOR 


Neutron stars are expected to have a solid inner crust 
of nonuniform neutron-rich matter above a liquid man¬ 
tle. The phase transition from solid to liquid is thought 
to be weakly first order and can be found by comparing 
a detailed model of the nonuniform crust to the liquid 
(see for example Q). Yet in practice, model calcula¬ 
tions yield very small density discontinuities at the tran¬ 
sition. Therefore a good approximation is to search for 
the density where the uniform liquid first becomes unsta¬ 
ble to small amplitude density oscillations (see for exam¬ 
ple H^]). This method would yield the exact transition 
density for a second order phase transition. 

The stability analysis of the uniform ground state is 
based on the relativistic random-phase-approximation 
(RPA) of Ref. |ilj for a system of electrons, protons, 
and neutrons. The approach is generalized here to ac¬ 
commodate the various nonlinear couplings among the 
meson fields. We start by considering a plane wave den¬ 
sity fluctuation of momentum q = |q| and zero energy 
go = 0. To describe small amplitude particle-hole (or 
particle-antiparticle) excitations of the fermions we com¬ 
pute the longitudinal polarization matrix that is defined 
as follows: 


Hr, = 


/ n oo 

0 

0 

0 

0 

H" + n§ 

n p 

m 

n n 

m 

0 

n m 

n p 

11 00 

0 

V 0 


0 

TT ra 

ii 00 


(4) 


Here the one-one entry describes electrons, the two- 
two entry protons plus neutrons interacting via scalar 
mesons, the three-three entry protons interacting with 
vector mesons and the four-four entry neutrons interact¬ 
ing with vector mesons. The individual polarization in¬ 
sertions are given by 


*n s (g, q 0 ) = 
Hl m (q,q 0 ) = 

zn 00 (g, go) = 


d A p 

W? 

d A p 

W 

d A p 


Tr 


Tr 


Tr 


G{p)G(p + g) 
G(p)joG(p + g) 
G{jp)"ioG(p + g)7o 


(5a) 

(5b) 

(5c) 


where Tr indicates a trace over Dirac indices. Note that 
the fermion Green’s function has been defined as 


G(p) = (0 + M*) 


1 


Z7T 


*2 _ + f *2 


+ — s(p* - K)e(k F - \ P \) 


( 6 ) 


Here k F is the Fermi momentum, M* = M — g s <t>o is 
the nucleon effective mass, E* = (p 2 + M* 2 ) 1 / 2 , and 
P*p = Pp~ (PvVpEgpbp/ 2) (with the plus sign for protons 
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and the minus sign for neutrons). Note that in the case 
of the electrons M* = m e and p p = p p . Explicit analytic 
formulas for IIoo, II S , and II m in the static limit (qo = 0) 
are given in the appendix. 

The lowest order meson propagator is computed in 
Ref. m in the absence of nonlinear meson couplings. It 
is given by, 


/ d 9 0 ~ d 9 0 \ 

0 -d° 0 0 

-d g 0 d g + d° v + d° p d°-d° p ] 

\0 0 d°-d° d° v + d° p / 


(7) 


Expressions for the photon and for the various meson 
propagators in the limit of no nonlinear meson couplings 
are given as follows: 


d a q 2 q 2 


C = 


9s 


q 2 + m 2 

,2 




ffv 


d" = 


q 2 + TO 2 

q 2 + m 2 


(8a) 

(8b) 

(8c) 

(8d) 


The appearance of a minus sign in the one-three element 
of D° l relative to the one-one element is because electrons 
and protons have opposite electric charges. 

The addition of nonlinear couplings in the Lagrangian 
leads to a modification of the meson masses. Effective 
meson masses are defined in terms of the quadratic fluc¬ 
tuations of the meson fields around their static, mean- 
field values (the linear fluctuations vanish by virtue of 
the mean-field equations). That is, 


TO 


*2 

s 


d 2 C * 2 _d 2 C 

dcj) 2 ' mv ~ + ^v 2, TOp 


d 2 C 

+ M' 


(9) 


This yields the following expressions for the effective me¬ 
son masses in terms of the static meson fields and the 
coupling constants defined in the interacting Lagrangian 
of Eq. (§): 

to* 2 = to 2 + g 2 ^/c$ 0 + - 2A s B$j , (10a) 

to; 2 = ml + gl gw 2 + 2A V R 2 ) , (10b) 

mf = m 2 p + g 2 p (2A S $ 2 + 2A V W 2 ) . (10c) 


Note that the following definitions have been introduced: 
$o = g s 4>o , W 0 =g v V 0 , and B 0 = g p b 0 . 

Further, the new couplings between isoscalar and 
isovector mesons (A s and A v ) lead to additional off di¬ 
agonal terms in the meson propagator. These arise be¬ 
cause the quadratic fluctuations around the static solu¬ 
tions generate terms of the form 


d 2 C 

dfodbo 


7^0 


and 


d 2 C 
dV 0 db 0 


^0. 


( 11 ) 


For simplicity we only consider here the following two 
cases: i) (A s ^ 0 and A v = 0) or ii) (A s = 0 and A v ^ 0), 
and neglect the (slightly) more complicated case in which 
both coupling constants are different from zero. 

For the first case of (A s ^ 0 and A v = 0) the new com¬ 
ponents of the longitudinal meson propagator become 


, = _ 9s (g 2 + ™* 2 ) _ 

(<f + m* 2 ){q 2 + to* 2 ) + (Ag s g p A s $ 0 B 0 ) 2 ’ 

, = _ (ffp/4)(g 2 + to. s * 2 ) _ 

P (q 2 + m* 2 ){q 2 + m* 2 ) + {4g s g p A s $ 0 B Q ) 2 ’ 

, __ 2g 2 g 2 A s $ 0 B 0 _ 

sp ( q 2 + m* 2 )(q 2 + m* 2 ) + (Ag s g p A s ® 0 B 0 ) 2 ' 


(12a) 

(12b) 

(12c) 

(12d) 


With these changes the modified longitudinal meson 
propagator now reads, 


D l 


dg 

0 

-dg 

0 

0 

d s 

dsp 

dsp 

-dg 

dsp 

dg + d v + d p 

d v — d t 

0 

dsp 

a 

-e 

1 

> 

dv + d t 


(13) 


Alternatively, for the case of (A s = 0 and A v ^ 0) we 
obtain, 


g 2 (q 2 + to.; 2 ) 

v ( q 2 + m* 2 )(q 2 + to; 2 ) - (Ag v g p A v W 0 B 0 ) 2 
= _ (ffp/4)(g 2 + to; 2 ) _ 

p (q 2 + m* 2 )(q 2 + to; 2 ) - (4g v g p A v W 0 B 0 ) 2 

, = _ - 2glg 2 A v W 0 B 0 _ 

vp ( q 2 + mf)(q 2 + m* 2 ) - {Ag v g p A v W 0 B 0 ) 2 


(14a) 
, (14b) 
, (14c) 
■ (14d) 


With these changes the modified longitudinal meson 
propagator now reads, 


D l = 


( d 9 

0 

-dg 
. 0 


0 

ds 

0 

0 


-dg 

0 

dg + d v + d p + 2d. 
d v d n 


vp 


0 

0 

dv d p 

r~\~d p 2d vp / 


(15) 


The uniform system becomes unstable to small ampli¬ 
tude density fluctuations of momentum transfer q when 
the following condition is satisfied: 


det 


1 - DL(q)RL(q,qo = 0) 


< 0 . 


(16) 
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We estimate the transition density (p c ) between the inner 
crust and the liquid interior as the largest density for 
which Eq. ( 0 ) has a solution. Our results for p c are 
listed in Tables 10 and also shown in Fig. [|. We find a 
strong correlation between the neutron skin of 208 Pb and 

□ 



FIG. 1. The transition density p c at which uniform matter 
becomes unstable to density oscillations as a function of the 
neutron skin in 208 Pb. The solid line is for the Z271 parameter 
set with A v ^ 0 while the dashed curve uses Z271 with A s ^ 0. 
The dotted curve is for the S271 set and the dot-dashed curve 
for NL3, both of these with A v ^ 0. 

At the lower densities of the inner crust the system is 
nonuniform and may have a very complex structure that 
may include spherical, cylindrical, and plate-like nuclei, 
bubbles, rods, plates, etc. [ fi 2 l|v} . At present we do not 
have microscopic calculations of these structures in our 
models. Therefore we adopt a simple interpolation for¬ 
mula to estimate the equation of state in the inner crust. 
That is, we assume a polytropic form for the EOS in 
which the pressure is approximately given by s 

P(e) = A + Be i/3 , (17) 

where e is the mass-energy density. The two constants 
A and B in Eq. 0 ) are chosen so that the pressure is 
continuous at the boundary between the inner crust and 
the liquid interior (determined from the RPA analysis) 
and at the boundary between the inner and the outer 
crusts. For the low density outer crust we assume the 
EOS of Baym, Pethick, and Sutherland(BPS) Q up to 
a baryon density of Pouter = 2.57 x 10 ~ 4 fm -3 which corre¬ 
sponds to an energy density of e 0 uter= 4.30 x 10 11 g/cm 3 
(or 0.24 MeV/fm 3 ) and a pressure of P ou ter = 4.87 x 
10 ” 4 MeV/fm 3 . Thus the two constants A and B of 
Eq. ( 0 ) are adjusted to reproduce P ou ter at e 0 uter and 
the pressure of the uniform liquid, calculated within the 
relativistic mean-field (RMF) approach, at e c which is 
the energy density corresponding to p c . That is, 

{ Pbps(c) , for ^min — € — ^outer 5 

A + Pe 4/3 , for Couter < e < e c ; (18) 

Prmf(c) , for e c < e . 


Note that e m i n = 5.86 x 10 -9 MeV/fm 3 is the minimum 
value of the energy density included in the equation of 
state. This value corresponds to a minimum pressure 
of Pmin = 6.08 x 10 ~ 15 MeV/fm 3 , which is the value 
at which we stop integrating the Tolman-Oppenheimer- 
Volkoff equations. That is, the radius R of a neutron 
star (see Sec. IV) is defined by the expression P(R) = 
Pmin- For the relativistic mean field interaction (TM 1 ) of 
Ref. |ji|, the relatively simple procedure presented here 
is a good approximation to the more complicated explicit 
calculation of the EOS in the inner crust Eel] . 


IV. RESULTS 

Figures [|-|| and Tables |l|-0 show the radii of neutron 
stars of mass 1/3,1/2,3/4 and 1.4 Mg as a function of 
the neutron skin ( R n -R p ) of 208 Pb. One might expect a 
strong correlation between the radius of a neutron star 
and the neutron radius of 208 Pb, as the same pressure of 
neutron rich matter that pushes neutrons out against sur¬ 
face tension in 208 Pb pushes neutrons out against gravity 
in a neutron star |2^] . However, the central density of a 
1.4 Mq neutron star is a few times larger than normal 
nuclear-matter saturation density po- Thus, P(1.4 Mq) 
depends on the EOS at low and high densities while R n 
only depends on the EOS for p<p 0 . Softening the EOS 
[i.e., decreasing the pressure) at high densities will de¬ 
crease P(1.4 Mq) without changing R n . Hence, while 
Fig. I shows a definite correlation—P(1.4 Mq) grows with 
increasing R n -R p — there is a strong model dependence. 

In contrast, the central density of a \Mq star is of 
the order of po, so P(|Mq) and R n depend on the EOS 
over a comparable density range. As a result, we find a 
strong correlation and weak model dependence in Fig. [^. 
For example, if R n -R p in 208 Pb is relatively large, e.g, 
R n —Rp ~ 0.25 fm, then P(^M 0 )« 16 km. Alternatively, 
if R n — R p w 0.15 fm, then R(^Mq) < 13 km. This is 
an important result. It suggest that even if observations 
of low mass neutron stars are not feasible, one could still 
infer their radii from a single nuclear measurement. Note 
that the results for a |M 0 neutron star (Fig. ||) follow a 
similar trend. 

We conclude this section with a comment on \Mq neu¬ 
tron stars. Parameter sets that generate very large val¬ 
ues for R n —R p have large pressures near po- This implies 
that the energy of neutron rich matter rises rapidly with 
density. In turn, this leads (because all parameter sets 
are constrained to have the same symmetry energy at 
p = 0.1 fm -3 ) to lower energies and lower pressures at 
very low density as compared with parameter sets with 
smaller values for R n — R p . This low-density region is 
important for low mass neutron stars and could explain 
why R(^Mq) in Fig. |2j actually decreases with increasing 
neutron skin for R n — R p > 0.23 fm. 
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FIG. 2. Radius of a neutron star of mass 1/3 Mg as a 
function of the neutron skin in 208 Pb. The solid line is for the 
Z271 parameter set with A v 4 0 while the dashed curve uses 
Z271 with A s 4 0. The dotted curve is for the S271 set and 
the dot-dashed curve for NL3, both of these with A v 4 0. 



FIG. 3. Radius of a neutron star of mass 1/2 Mg as a 
function of the neutron skin in 208 Pb. The solid line is for the 
Z271 parameter set with A v 4 0 while the dashed curve uses 
Z271 with A s 4 0- The dotted curve is for the S271 set and 
the dot-dashed curve for NL3, both of these with A v 4 0. 


FIG. 4. Radius of a neutron star of mass 3/4 Mq as a 
function of the neutron skin in 208 Pb. The solid line is for the 
Z271 parameter set with A v 4 0 while the dashed curve uses 
Z271 with As 4 0- The dotted curve is for the S271 set and 

t/hp rlnt.-rlashprl nirvp fnr NT,3 Tnt/h r»f t/hpsp wit/h A.. zA 0 



FIG. 5. Radius of a neutron star of mass 1.4 M© as a 
function of the neutron skin in 208 Pb. The solid line is for 
the Z271 parameter set with A v 4 0 while the dashed curve 
uses Z271 with A s 4 0- The dotted curve is for the S271 set 
and the dot-dashed curve for NL3, both of these with A v 4 0. 
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TABLE II. Results for the NL3 parameter set with A s = 0. The NNp coupling constant g 2 and the neutron minus proton 
root mean square radius for 208 Pb (in fm) are given along with the radii of 1/3, 1/2, 3/4 and 1.4 Mq neutron stars in km. 
Finally, the transition density p c between the inner crust and liquid interior is given in fm~ 3 . 


A v 

9p 

R n -R P ( 208 Pb) 

R(iM 0 ) 

R(|M 0 ) 

R(| M q ) 

R(1.4M 0 ) 

Pc 

0.030 

127.0 

0.1952 

16.766 

14.789 

14.142 

14.175 

0.0854 

0.025 

115.6 

0.209 

18.37 

15.59 

14.60 

14.38 

0.0808 

0.020 

106.0 

0.223 

19.49 

16.15 

14.93 

14.52 

0.0746 

0.015 

97.9 

0.237 

19.73 

16.39 

15.10 

14.61 

0.0675 

0.010 

90.9 

0.251 

19.31 

16.40 

15.20 

14.68 

0.0610 

0.005 

84.9 

0.265 

18.70 

16.35 

15.31 

14.81 

0.0558 

0.000 

79.6 

0.280 

18.10 

16.27 

15.47 

15.05 

0.0519 


TABLE III. Results for the S271 parameter set with A s = 0. The NNp coupling constant g 2 p and the neutron minus proton 
root mean square radius for 208 Pb (in fm) are given along with the radii of 1/3, 1/2, 3/4 and 1.4 Mq neutron stars in km. 
Finally, the transition density p c between the inner crust and liquid interior is given in fm~ 3 . 

A v 

Sp 

R n - R P ( 208 Pb) 

R(|M 0 ) 

R(|m q ) 

R(|M 0 ) 

R(1.4M 0 ) 

Pc 

0.05 

127.8389 

0.1736 

15.88 

14.06 

13.43 

13.25 

0.0923 

0.04 

116.2950 

0.1895 

17.75 

15.00 

13.96 

13.47 

0.0866 

0.03 

106.6635 

0.2054 

19.25 

15.77 

14.42 

13.65 

0.0794 

0.02 

98.5051 

0.2215 

19.80 

16.24 

14.76 

13.82 

0.0717 

0.01 

91.5061 

0.2378 

19.54 

16.46 

15.08 

14.07 

0.0648 

0.00 

85.4357 

0.2543 

18.95 

16.53 

15.43 

14.56 

0.0594 


TABLE IV. Results for the Z271 parameter set with A s = 0. The NNp coupling constant g 2 and the neutron minus proton 
root mean square radius for 208 Pb (in fm) are given along with the radii of 1/3, 1/2, 3/4 and 1.4 Mq neutron stars in km. 
Finally, the transition density p c between the inner crust and liquid interior is given in fm~ 3 . 

A v 

9p 

Rn 

- R p ( 208 Pb) 

R(|M 0 ) 

R [\Mq) 

R(fM 0 ) 

R(1.4M 0 ) 

Pc 

0.14 

139.3368 


0.1525 

13.799 

12.709 

12.293 

11.616 

0.0974 

0.12 

129.2795 


0.1650 

15.012 

13.379 

12.688 

11.748 

0.0936 

0.10 

119.5245 


0.1771 

16.219 

14.039 

13.080 

11.880 

0.0890 

0.08 

112.9710 


0.1900 

17.405 

14.696 

13.481 

12.016 

0.0845 

0.06 

106.2682 


0.2026 

18.31 

15.28 

13.89 

12.18 

0.0796 

0.05 

103.2065 


0.2090 

18.61 

15.53 

14.09 

12.29 

0.0771 

0.04 

100.3162 


0.2154 

18.82 

15.76 

14.30 

12.43 

0.0747 

0.03 

97.5834 


0.2218 

18.95 

15.96 

14.53 

12.62 

0.0725 

0.02 

94.9956 


0.2282 

18.98 

16.12 

14.75 

12.89 

0.0703 

0.01 

92.5415 


0.2347 

18.94 

16.25 

14.98 

13.27 

0.0683 

0.00 

90.2110 


0.2413 

18.88 

16.36 

15.20 

13.77 

0.0665 
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TABLE V. Results for the Z271 parameter set with A v = 0. The NNp coupling constant g 2 and the neutron minus proton 
root mean square radius for 208 Pb (in fm) are given along with the radii of 1/3, 1/2, 3/4 and 1.4 M 0 neutron stars in km. 
Finally, the transition density p c between the inner crust and liquid interior is given in fm -3 . 


A s 

9 2 p 

An - A P ( 208 Pb) 

R(iM 0 ) 

R(|M 0 ) 

R(|M 0 ) 

R(1.4M 0 ) 

Pc 

0.06 

146.6988 

0.1640 

13.34 

12.71 

12.52 

11.98 

0.0844 

0.05 

132.8358 

0.1775 

14.63 

13.48 

13.01 

12.19 

0.0830 

0.04 

121.3666 

0.1907 

15.99 

14.27 

13.51 

12.42 

0.0807 

0.03 

111.7205 

0.2036 

17.28 

15.02 

14.00 

12.68 

0.0777 

0.02 

103.4949 

0.2163 

18.27 

15.66 

14.46 

12.97 

0.0741 

0.01 

96.3974 

0.2288 

18.79 

16.10 

14.86 

13.32 

0.0702 

0.00 

90.2110 

0.2413 

18.88 

16.36 

15.20 

13.77 

0.0665 


V. DISCUSSION AND CONCLUSIONS 

A number of relativistic effective field theory parame¬ 
ter sets have been used to study correlations between the 
radii of neutron stars and the neutron radius R n of 208 Pb. 
An RPA stability analysis was employed to find the tran¬ 
sition density between the nonuniform inner crust and 
the uniform liquid interior. For the nonuniform outer 
crust we invoked the EOS of Baym, Pethick, and Suther¬ 
land 0. Then, a simple polytropic formula for the EOS, 
approximately valid for most of the crust 0], was used 
to interpolate between the outer crust and the liquid in¬ 
terior. This simple, yet fairly accurate, procedure allows 
us to study the EOS for a variety of parameter sets that 
predict a wide range of values for the neutron radius of 

2°8p b . 

For a “normal” 1.4M 0 neutron star we find central 
densities of several times normal nuclear matter satura¬ 
tion density (^0^0.15 fm 3 ). Because the neutron radius 
of 208 Pb does not constrain the high-density component 
of the EOS, we find a strong model dependence between 
the radius of a 1.4M 0 neutron star and R n . In contrast, 
the central density of a low mass neutron star is close to 
Pq. Therefore properties of the low mass star are sensi¬ 
tive to the EOS over the same density range as R n . As 
a result, we find a strong correlation and a weak model 
dependence between the radius of a 0.5M 0 neutron star 
and R n (see Fig. |j). Thus, it should be possible to infer 
some properties of low mass neutron stars from detailed 
measurements in atomic nuclei. 

Understanding the density dependence of the equation 
of state is particularly interesting. A softening of the 
EOS at high density (where the pressure rises slower than 
expected) could signal the transition to an exotic phase, 
such as pion/kaon condensates, strange quark matter, 
and/or a color superconductor. Yet obtaining definitive 
results on the density dependence of the EOS may require 
measuring the radius of neutron stars for a broad range 
of masses. This may be difficult as most compilations to 
date find neutron star masses in the very narrow range of 
1.25 — 1.44 Mq (t|. Further, important ambiguities will 


remain if the radius of a 1.4M 0 neutron star proves to be 
moderately small. Would a small radius be an indication 
that the EOS is relatively soft at all densities and there 
is no phase transition, or is the EOS stiff at low density 
and undergoes an abrupt softening at high density from 
a phase transition? 

Therefore, it is important to also make measurements 
that are exclusively sensitive to the low density EOS. One 
obvious possibility is the radius of a low mass (« O.5M 0 ) 
neutron star as it central density, close to po, is much 
smaller than that of a 1.4M 0 star. However such low 
mass stars may be very rare because they are hard to 
form. Notably, the neutron radius of a heavy nucleus, 
such as 208 Pb, contains similar information. Indeed, we 
find a strong correlation and a weak model dependence 
between the neutron radius of 208 Pb and the radius of a 
0.5 Mq neutron star. This allows one to use nuclear in¬ 
formation to infer the radius of a low mass neutron star. 
Hence, comparing this inferred radius to the measured ra¬ 
dius of a ~ 1.4M 0 neutron star, should provide the most 
complete information to date on the density dependence 
of the equation of state. 

This work was supported in part by DOE grants DE- 
FG02-87ER40365 and DE-FG05-92ER40750. 






VI. APPENDIX 


The polarizations are defined in Eqs. ( |5a| , |5b|j5c[ ) and describe particle-hole or particle antiparticle excitations. In the 
static limit (energy transfer go = 0) the scalar polarization for momentum transfer q is given by 


n.(5,o) = s? 


kpEjr - 


3M * 2 + ^)\n kF + EF ' 2EfE 


M* 


In 


2 kp — q 

2E 3 4 * * * * * * * 12 13 14 15 16 

qEp — 2 kpE 

‘Ikj? q 

9 

qEp + 2 kpE 


(19) 


with Fermi momentum kp, nucleon effective mass M*, Ep = (kp + M* 2 ) 1 / 2 , and E = (q 2 /4 + M * 2 ) 1 / 2 . Likewise, 
the longitudinal polarization is given by 


IIoo(9, 0) — —- 


hpEp-Un^±^-^ 

36 M* 3 q 


(M 


*2 


kp 


3 2N 1 
-4 9 )ln 


2kp — q 


2 kp + q 


!< M' 2 -y) h 


qEp — 2 kpE 


qEp 


2 kpE 

( 20 ) 


while the mixed scalar-vector polarization becomes 


M* 




n m(g, °) = ^ <j kp - (— - 7 ) In 


2kp — q 


2 kp + q 


( 21 ) 
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